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We derive relations among branching fractions in the exclusive decay of charmonia to light flavour 
meson pairs assuming factorization between the quark spin and spatial degrees of freedom. With the 
further assumption that these amplitudes can be described by flux-tube models, we assess prospects 
for production of hybrid mesons in charmonium decays. 

I. INTRODUCTION 

o 

i^« ■ The decays of tp to light hadrons via qq intermediate OZI diagrams have been investigated by several authors 

[1, 2]. These analyses showed that while leading OZI diagrams dominate and give a good first description of data, 
more detailed fits implied that disconnected diagrams cannot be neglected in general[2]. These studies primarily 
related states of different flavour but having the same spin and spatial quantum numbers. To go beyond that it is 
neccessary to make an assumption concerning the intermediate qq state. Even within the leading OZI assumption, 
£f) ' such efforts have foundered on the need to model the radial wavefunction of the intermediate qq state; the lack of 
£N| , knowledge of which radial excitations ('«') dominate, and the need to model the specific wavefunctions, introduced 
model dependences that severely limited useful predictions for the final states. 

Recently we have shown [3] that results from lattice QCD imply a factorization of L and S in strong decays decays, 

O . and that the created pair are spin triplet. This hypothesis factors out the spatial degrees of freedom and gives relations 

( "L-S" relations) among amplitudes for the decays of states sharing the same spatial quantum numbers but differing 

jr~^, in spin and total angular momentum that are independent of of the spatial wavefunctions. In this paper we apply 

,-h such relations to the hadronic decays of charmonia to light flavour meson pairs, assuming that such decays can be 

modelled as the strong decay of an intermediate light quark qq state. 

By exploiting the factorization of amplitudes in OZI strong decays, we can avoid the n-dependence problem in so 
far as we leave the spatial part of the decay amplitude undetermined, and instead exploit the relations between decay 
amplitudes arising from the factorisation of L and S. Such relations apply in a symmetry limit, provided there is not 
yy-Z a double conspiracy of a single n dominating and kinematic node for that n. Within this more modest assumption 
£ — . we show that if the transition is indeed driven by a virtual qq state, as implicitly assumed in analyses [1, 2], then 
C*~) ' some qualitative and semi-quantitative conclusions can be drawn about its nature and new relations among some final 
states obtained. 

We find that data rule out an intermediate qq hybrid state and appear to be consistent with a coupling of cc-^ 
f*"*. \ qq where the qq is 'canonical', in the sense of having the same 2S+1 Lj configuration as the initial cc state. If 
^^> confirmed by further data this will provide a challenge for dynamical models to explain. It suggests that hybrid 
qq production at short distances is suppressed, which would imply that their production from electromagnetic or 
weak currents, e.g. in e + e~ or B decays will be feeble. 

This leaves open the question of whether hybrids can be produced in charmonia decays in association with other 
hadrons. Within our model there is a selection rule suppressing the production of a negative parity hybrid along with 
a IS meson; this applies to certain modes which a priori looked promising, such as \i — + 7T7Ti . On the other hand, we 
find that other exotic hybrids such as the 2^ 62 may be produced in sufficient measure to be observed, for instance 
in tp — ► &27T- 

In section II we set up the model. In section III we confront the model predictions with data on charmonia decays. 
In section IV wc apply the model to the production of hybrid mesons. 
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II. THE MODEL 

We are interested in the decay topology in which the cc annihilate into gluons or a virtual photon, creating a light 
quark qq pair which then decays by the creation of an additional light quark pair QQ, 

cc^ qq — > qQ + Qq, (1) 

as depicted in FIG. 1. Flavour disconnected decays cc — > qq — > qq + QQ are not negligible but are perturbations on 
the main phenomenology and will be ignored for this first survey. 
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FIG. 1: cc annihilation into a virtual qq state which decays by the creation of a QQ pair. 

We postulate that the intermediate qq state has well defined j = I © s and can be described as a non-relativistic 
quark- flux tube system, a meson m qq , which decays by flux tube breaking in the usual way to mesons ml + m^. 
Thus we describe the transition 

m cc ^ m qg ^ m iQ + m Qv ( 2 ) 

by the amplitude 

(mfm$ q \V\m d5 } = ^ (mfm$ q \H I \mf)(mf\V Q cD\m c ' 5 ) (3) 

\mf)(mf\ 

where 

• l m ? <? )( rn f 9 | i s a complete set of quark-flux tube states m qq consistent with the J PC quantum numbers, 

• (tti^ \Vqcd\it^ cc ) is the amplitude for creation of the virtual state m qq via cc — ► ggg("f*) — * qq, and 

• {m\ m^ \Hi\ m i) is the amplitude for the decay of the virtual state m gq via qq — > qQ + Qq. 

Writing eqn. (3) in another way, the full decay amplitude is proportional to the strong decay amplitude of an 
intermediate state \m qq ), 

{mfmf ! \V\m cZ ) = {mfm^H^m™) (4) 

where \m qq ) is a linear combination of possible configurations m qq with mixing angles governed by the Vqcd interac- 
tion, 

K 3 >= E \mf){mf\V QCD \m c - c ) (5) 

\mf)(mf | 

The complete set of intermediate virtual qq states \m qq )(i r n qq \ includes all possible orbital, radial and gluonic excita- 
tions consistent with the initial state quantum numbers, and a priori the hierarchy of mixing angles {m qq \VQCD\ i nn, cc ) 
between such states is unknown. The amplitude (w.* 9 |VQC£)|m cc ) describes a complicated non-pcrturbative pro- 
cess which we do not attempt to calculate; instead, we suppose that there is a unique choice of m qq for which 



(m qq i\VQCD\m cc ) is large; that is to say, the wavefunction of \m qq ) is dominated by a qq state \m qq ) in a unique 

S,L and A configuration. With this assumption, the strong decay matrix element (mf rri2 \Hi\m qq ) relates the 
decay amplitudes of different final states using eqn. (4), subject only to the factorisation of spin and orbital angular 
momentum quantum numbers. 

We will use the notation of Ref [4] . A quark- flux tube system can be described by a ket of the form 

\{s®nAl p )j) (6) 

where s, I and j are the usual spin, orbital and total angular momentum quantum numbers, n is the quark radial 
quantum number, A is the gluonic angular momentum and p is a parity label appropriate to A ^ states which 
occur in parity doublets. In anticipation of our proposal that the intermediate virtual state is dominated by a single 
configuration, let us associate \m qq ) with a unique ket of the above form 

\m qq ) = \{s®nAP) 3 ). (7) 

The final state (m\ wij | consists of a pair of mesons likewise defined, coupled to j\i and in a relative partial wave 
L, in turn coupled to angular momentum j, 

(mfm^\ = ((((«! <g) mAiif ) h <g> (s 2 ® n 2 A 2 l% 2 ) J2 ) jl2 <g> L) jm \. (8) 

In the flux tube model the decay operator Hi is given by 

Hi = xi ■ V (9) 

where xi creates a spin 1 QQ pair and V acts on the quark degrees of freedom. Thus we are interested in matrix 
elements of the form 

(mfm^lHjlm™) = <((( Sl <g> mAil? 1 ) h <8> {s 2 <8> ra 2 A 2 Z£ a ) i2 ) jia ® L^Hxi • V||(s <8> nAP),). (10) 

The decay rate is proportional to the sum over the squared transition amplitudes for all possible couplings j\ 2 and 
partial waves L; it is useful to define a corresponding bracket 

[mfm^\m qq ] = ^ ((((si ® mAiZf ) J - 1 ® (s 2 ® ^ 2 A 2 /f 2 ) j2 ), 12 ® L),-|| X i • V||(* ® nAP),-) 2 . (11) 

Following eqn. (4), the branching fraction for a charmonium state to decay to a light meson pair ml + m^ is 

6.r.(?7i" -» mp + mf 1 ) oc pSiSatmf^m^lm 45 ] (12) 

where E\ , i? 2 are the energies of mesons with mass M\ , M 2 produced with momentum p. 

Our starting point will be the derivation of the matrix element (10) as presented in refs. [4, 5]. In the approach 
presented there, the first step is to separate the spin and space degrees of freedom. This involves a recoupling of the 
states of good ji,j 2 ,ji 2 into states of good s\ 2 , l\ 2 , If, whence the spin matrix element is simply given in terms of a 
9-j coefficient. In this way the full matrix element 

((((s 1 ®niA 1 % 1 ) jl ®(s 2 ®n 2 A 2 l 2 ,2 ) h ) jl2 ®L) j \\xi-V\\(s®nAl p ) j ), (13) 

is expressed as linear combination of spatial matrix elements of the form 

(((niAiJ? 1 ®n 2 A 2 l p 2 2 ) ll2 ®L) lf \\V\\nAP). (14) 

This leads to relationships among amplitudes to different final states sharing the same spatial quantum numbers but 
differing in spin and total angular momentum. Such relations apply in the limit that the final states under comparison 
have the same masses, radial wavefunctions and decay momenta. The relations are also independent of the radial 
wavefunction of the initial state, hence they translate directly into relations among the branching fractions using eqn. 
(12) provided that there is not a double conspiracy in which m qq is dominated by a single radial configuration n which 
has a kinematic node suppressing certain channels. 
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III. L-S RELATIONS AMONG AMPLITUDES 

We now turn to the question of the intermediate state \m qq ), which can in general be a linear combination of states 
\m? 9 ) consistent with the J PC quantum numbers of \m cc ), as in eqn (5). For the charmonia states of interest, the 
angular momentum configurations consistent with the J quantum numbers are [5] : 

(15) 
(16) 
(17) 
(18) 
(19) 

where in the above we have used the parlance of molecular physics, in which the quantum numbers I = 0, 1, 2, . . . are 
labelled S,P,D. . . and A = 0, 1, 2 . . . arc analogously labelled E, II, A . . .. 

In the following sections we discriminate among the possible angular momentum configurations s, A and / for the 
intermediate states m qq with reference to data on rj c , x an d ip decays. In ref. [5] it was shown that the selection rule 
of refs. [6, 7] forbidding the decay of IIP hybrids to identical S-wave mesons applies also to IID and AD-type hybrids. 
Using this, we can rule out hybrid configurations in r\ c and x decays which show prominent production of identical 
S-wave pairs. In the case of ip the spin singlet selection rule argues against the hybrid intermediate state. 

We hypothesise that m qq is dominated by the E-state with the smallest I consistent with the quantum numbers; 
we refer to this as the "canonical choice" wavefunction. Within this hypothesis we express the production rates as 
linear combinations of spatial matrix elements of the type (14) with the superfluous A and p labels dropped, 

(((ni/i <gi n 2 l 2 )i 12 ®L) lf \\V\\nl). (20) 

The resulting L-S relations arc consistent with data and make predictions for modes which have not yet been observed. 

A. Tj c decays 

The virtual qq state m qq corresponding to the r\ c could have ^^ESo or 3 IIPo assignments. On heuristic grounds we 
expect the latter to be suppressed owing to a node in both the quark and flux tube wavefunctions at the origin. If the 
virtual qq corresponding to the r\ c were dominantly 3 ITPo, the production of identical 3 Si pairs would be forbidden 
not only by the standard hybrid selection rule forbidden the decay to identical S-wave pairs [6, 7] but also by a more 
general selection rule derived in ref. [5], 

[ 3 Si 3 Si| 3 inv] = o, (21) 

which relies only on L-S factorisation. In this context it is notable that three vector-vector modes have been observed 
with considerable branching fractions [8], 

b.r.{r) c -> pp) = (20 ± 7) x 10~ 3 (22) 

b.r.(r} c ^K*K*) = (9.2 ± 3.4) x 10~ 3 (23) 

b.r.(r) a -»• #) = (2.7 ±0.9) x 10~ 3 (24) 

These branching fractions are comparable with multi-body modes 

b.r.(r] c -> t/tttt) = (41 ± 7) x 10" 3 (25) 

b.r.^c -> K*°K-n + + c.c.) = (20 ± 7) x 10~ 3 (26) 

b.r.(r) c -> K* Q ~K*°ir + Tr-) = (15±8)xl0~ 3 (27) 

b.r.(r] c -> <j)K + K-) = (2.9 ± 1.4) x 10~ 3 (28) 

and for pp, at least, exceed those of the 1P+1S modes expected to dominate for decay of 3 IIPo state, 

b.r.(j] c -»■ a 2 7r) < 20 x 10~ 3 (29) 

b-r.(Vc -» M) < 11 x 10" 3 (30) 
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This supports the conjecture that m qq will be dominated by the ^"ESo rather than 3 IlPo configuration. Assuming 
therefore the canonical choice '"ESo wavcfunction dominates, one can express the allowed decay amplitudes in terms 
of common spatial matrix elements and therefore extract relations between them. For decays to S-wave pairs, only 
vector- vector pairs are allowed. For P+S final states there are several possible decays, each of which can be expressed 
uniquely in terms of S- and D-wave spatial matrix elements Si and D\ defined 

Si = (((niP®naS)i®S)i||V||nS) (31) 

£>i = (((niP®naS)i®D)i||V||nS) (32) 

In the above we leave the radial quantum numbers n, n\ and ni explicit as a reminder of the fact that the relations we 
derive are independent of the radial wavcfunctions of the initial and final states. If specific models are subsequently 
developed, the resulting wavcfunctions can be applied to these immediately. The appropriate decompositions are 
derived in ref. [4] and are presented explicitly in ref. [5]. For the width bracket of the form (11), the result is 

(33) 

(34) 

(35) 
and thus 

[ 3 Po 1 SoI 1 S ] + [ 3 P 2 'SofSo] = [ x Pi "SifSo] (36) 

Thus the decay rates for isovector modes are related 

[a 7r|?7 c ] + [a 2 7r|?7 c ] = [bip\r) c ] (37) 

At present there is an upper limit on the a27r mode [8]; if this mode can be isolated then eqn. (37) relates the 
branching fractions of b\p and a^it and this prediction could be confronted with experiment. For isoscalar pairs the 
analogous relationships follow immediately from the above with appropriate mixing angles. For mixed K±-K[ states 

\K X ) = cos^EPO+sin^EPi) (38) 

\K[) = -sin^EP^ + cos^EPi). (39) 

only the spin singlet part contributes to the decay amplitude, hence we predict 

[K K\r) c ] + [K 2 K\r] c } = [K^K*^] + [K^K*^]. (40) 

B. xo decays 

For xo decays, the only quantum number assignment available for m qq is 3 EPo- Proceeding as before, we use the 
decompositions presented in ref. [5] to express the allowed decay modes in terms of common spatial matrix elements. 
Both pscudoscalar and vector pairs are allowed, and the corresponding amplitudes are expressed in terms of two 
spatial matrix elements S and D, defined 

S = (((niS®n 2 S)o®S)o||V||nP) (41) 

D = (((niS®n2S)o®D) 2 ||V||nP) (42) 



The decompositions are 



from which we obtain the constraint 



pSo ^ofPo] = ^ (43) 

[ 3 Sl 3 Sl | 3 P ] = ^p- (44) 



[ 3 Si 3 S!| 3 P ] > l/S^So 'Sol'Po] (45) 



There are large phase space and momenta differences for 3 Si 3 Si and 1 So x So modes, and to quantify these would 
require specific model wavef unctions. The experimental modes [8] 



b.r.(xo -> K*°K*°) = 1.8 ± 0.6 x 1CT 3 (46) 

b.r.(xo^ K+R-) = 5.4±0.6xl(T 3 (47) 

and also the wwjryry and mr are nonetheless consistent with the above. 

Decays to P+S final states can be expressed in terms of a spatial matrix clement Pi defined as 



and the decompositions are 



Thus we predict 



and likewise for strange pairs 



P 1 = (((niP®n 2 S) 1 ®P)i||V||nP) (48) 



[ 3 P! ^ol'Po] = fPi 3 Si| 3 Po] = ^-- (49) 



[ai7r|Xo] = [hp\xo] (50) 



[KiKfoo] + [KiK\ Xo ] = [KiK*\ X o] + [K[K*\ X o] (51) 

These will be a useful experimental test of our hypothesis. 

C. xi decays 

For the 1 ++ sector, the m 9 ' assignments could be the canonical 3 XPi, or hybrid configurations -"-IIPi, 3 ADi. The 
first two have the same penalty for the quark wavefunction at the origin, but it is anticipated that the 3 SPi will win 
over the 1 IIPi because the latter will also have a penalty for the creation of a flux tube at the origin: it is not possible 
to have a transverse excitation in a flux tube with zero size. 

The only two-body mode that has been directly observed is K*K*: that this has comparable strength to multi-body 
modes argues against the 1 IlPi, 3 ADi assignments since the production of these identical particles would be forbidden 
owing to the (generalised) hybrid selection rule of ref . [5] . Hence the hybrid intermediate states are ruled out and we 
consider the L-S relations among P+S modes assuming the canonical 3 EPi dominates. The allowed decay amplitudes 
are expressed in terms of the spatial matrix clement P\ defined in cqn. (48) and three further spatial matrix elements 
P , P 2 and F a , 

P = (((niP<gm 2 S)i®P)o||V||nP) (52) 

P 2 = (((niP®n 2 S)i®P) 2 ||V||nP) (53) 

F 2 = (((n 1 P®n 2 S)i®F) 2 ||V||nP) (54) 

The decompositions are 

[3 p o i So|3pi] = 4P 2 + 3P 1 2 + 2x/I5P 1 P 2 + 5P 2 2 -4Pq(V3P 1 + x/5P 2 ) 

flBlfl „ D1 16Po 2 + 3P 1 2 -2yi5P 1 P 2 + 5P 2 2 +P (-8V3P!+8V5P 2 ) 

[ f 1 bo | fi\ - — (ob) 

l3n lc l3n , 80P 2 + 15P 1 2 -2\/l5P 1 P 2 + P 2 2 + 8Po (5V3P!-V5P 2 )+36P 2 2 

I p 2 5>o| PiJ - — (57) 

[lpi 3 Si| 3 Pi] = 4Po 2 + Pi 2 + P2 2 +F 2 * (58) 

from which we obtain the relation 

pPi 3 Sx | 3 Px] = [ 3 P ^ol'Pi] + [ 3 Pi ^ol'Pi] + [ 3 P 2 'SopPi] (59) 



This implies 

[ph\xi] = [tooIxi] + [ttoiIxi] + K°2|Xi] ( 60 ) 

and 

[ITIdlxi] + [K*K[\ X1 ] = [KKo\xi] + [KKtlxi] + [KK[\ X i} + [KK 2 \ X i] (61) 

Data on \i decays are sparse and we urge that tests of the above relations be investigated. 

D. X'2 decays 

For X 2 decays there are many possible assignments for the virtual qq state: two conventional states 3 £P2, 3 SF2; 
and three hybrids 1 ±1D2, 3 AD2, 3 AF2. The empirical prevalence of [8] 

X2 -* (/)(/), UJW,TTTT, r\T) (62) 

once again rules out hybrid interpretations. The two £ configurations can be discriminated by considering the L-S 
relations between pscudoscalar and vector pair amplitudes. For 3 P2 the decays can be expressed in terms of the 
spatial matrix elements S and D of eqns. (41) and (42), 

[% 1 S | 3 P 2 ] = ^ (63) 

[3 Sl 3 Sl | 3 P 2 ] = bS2 + 2D2 (64) 
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and hence there is a constraint 



By contrast, for 3 F2 one obtains 



[ 3 S 1 3 S 1 | 3 P 2 ]>-[ 1 So 1 So| 3 P 2 ]. (65) 



[ 3 S 1 3 S 1 | 3 F 2 ] = ^-[ 1 So 1 S | 3 F 2 ]. (66) 



Experimentally [8] 



b.r.( X 2 -> K* a ~K^) = (38 ± 9) x 1(T 4 (67) 

b.r.( X 2 -> K+R-) = (7.7 ±1.4) x 1CT 4 (68) 

b.r.( X2 -> K 8 K a ) = (6.7 ± 1.1) x 1(T 4 (69) 

which rules out the 3 F2 and once again the canonical configuration, in this case 3 P2, dominates. For isovectors, if wc 
associate the 2(7r + 7r~) mode with pp, the experimental result concurs: 

b-r.(x2 — ► PP) ~ 6 x b.r.(x2 — > Kit) (70) 

Amplitudes for the P+S modes are defined in terms of the spatial matrix elements Pi, P-i and F 2 of eqns. (48), (53) 
and (54), and the decompositions are 

[■P, 3 Sl|3 p 2] . »*' + .ff' + *') (73) 

This yields a relation 

[ 3 Pi 'So^Pa] + [ 3 P 2 1 S | 3 P 2 ] = fPi 3 Sx| 3 P 2 ] (74) 
which can be tested experimentally 



E. ip decays 

For ip the possible intermediate states mf 1 are 3 £Si, 3 £Di, 1 lTPi, 3 ADi. That the qq pair is created at a point 
suggests 3 ESi will dominate, and this appears to be confirmed by the data: any significant II or A admixture would 
be in conflict with known prevalence of experimental modes such as 

e + e~ — ► T] c ip, and (75) 

ip — > pir (76) 

owing to the standard selection rule, although this rule may be broken quite significantly for final states with different 
spatial wavcfunctions [9]. A stronger rule appears assuming only the factorisation of L and S in the form of the 
spin-singlet selection rule, which forbids decays of the type 

spin — ► spin + spin 0. (77) 

Experimentally, the decay 

V-^&lTT (78) 

is large; this would be forbidden if the virtual qq is in the hybrid IIPi state. 

With the canonical choice of 3 Si dominance we can express the decay rates to P+S states in terms of the spatial 
matrix elements Si and D\ of eqns. (31) and (32), 

I 1 ?! ^ol'Sx] = ^^ (79) 

fPo'Sxl'Si] = ^- (80) 

[ 3 Pi 'S^Sx] = 4Sl \ + 2 Dl2 (81) 

[ 3 P 2 3 S 1 | 3 S 1 ] = ^- (82) 

(83) 



and hence there are two independent relations 



pPi ^o^Si] = i[ 3 P 3 Si| 3 Sx] + i[ 3 P 2 "SxI'Sx] (84) 
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[ 3 Pi 3 Si| 3 Si] = ^Po'Si^Sil + ipPa^^Sx] (85) 

Thus for tp decays to isovector pairs the model predicts 

[MV>] = \\^p\^\ + \\aM^] (86) 

[aiPM = g[aopM + g[o2/<#] (87) 

Experimentally the modes a 2 p and b\-n have been measured 

b.r.{ip -> bin) = (5.3 ± 0.8) x 1CT 3 (88) 

b.r.(tp -+ a 2 p) = (10.9 ± 2.2) x 10~ 3 (89) 

and are in agreement with the weak constraint from the above 

MM > \[a*p\il>\> (90) 

6 

Using eqns. (86) and (87) gives a prediction for the scale of the as yet unseen a^p and a\p modes. As a first estimate 
we approximate the experimental result by [£>i7r|^0] = ^[a 2 |p] and thus from eqns (86) and (87) 

[aop\ip] = 2[6 l7 r|V] (91) 

[aip\il>] = 3[&i7#] (92) 



hence we predict 

b.r.(iP^a p) w lOxlO -3 (93) 

b.r.(ip -> oip) « 15xl0~ 3 (94) 

These modes will be difficult to identify but feed multi-body channels which are observed to be large, 

a\p — ► (pir)p — > 57r (95) 

ao/5 — ► (ujTnr)p — > w47r (96) 

The W7r + 7r + 7r~7r~ mode is observed with branching fraction (8.5 ± 3.4) x 10 -3 which is consistent with the above 
prediction. 

For strange pairs are the analogues of eqns (86) and (87) are 

[K^H'] + [K[K\i>] = ±[K K*\tP] + ^[K 2 K*\tP] (97) 

[KiK*\4] + [Kiir\1>] = ^[K K*\^] + ^[K 2 K*\i/>] (98) 

Experimentally 

b.r.(ip -> K*°K$ + c.c.) = (6.7 ± 2.6) x 1CT 3 (99) 

&.r.(V>-> JsTi(f400) ± A' T ) = (3.8 ±1.4) x 10~ 3 (100) 

&.r.(^-> A' 1 (1270) ± A' T ) < 3.0 x 10" 3 (101) 

If the Ai(1270) ± A" =F mode can be measured then the KqK* branching fraction is predicted by eqn. (97). 

Analysis of the 3 SDi possibility is in general model dependent. There is however a potential discriminator in the 
helicity selection rule of ref[3], which states that for 3 SSi the amplitude for decay to vector and tensor qq , such as 
tuf 2l vanishes when the tensor has helicity ±2. For 3 EDi by contrast, this amplitude is not suppressed. 

IV. HYBRID MESON PRODUCTION 

In this section we consider the prospects for hybrid production in charmonia decay. We take as our starting point 
the hypothesis confirmed in the previous section that the intermediate state m qq is dominated by a conventional £ 
state. With this assumption we have expressed decay rates to conventional meson final states with spatial quantum 
numbers ri\l\ and n 2 l 2 as linear combinations of spatial amplitudes of the form 

((Mi®n 2 Z 2 )i 12 ®£);/i|V|H. (102) 

Since these amplitudes are in general independent it is not possible to relate branching fractions among states with 
different spatial quantum numbers n\l\ and 712/2- 

Following the same approach, the decay amplitude to final states with spatial quantum numbers niAi^ 1 and 
"2A2/2 2 can be expressed in terms of spatial amplitudes of the form 

{{{n^ll 1 ®n 2 K 2 l p 2 2 ) h2 ® L) lf \\V\\nl). (103) 

The novel feature of the approach presented in refs. [4, 5] is that all spatial amplitudes of the form (103) are linearly 
related to spatial amplitudes of the form (102), leading in some cases to direct relationships between hybrid production 
amplitudes and those of conventional mesons. 

With this approach, it is possible to make general statements about hybrid production rates without having to 
make further assumption concerning the virtual state m q ^ . We restrict the discussion here to the production of a 
HIP hybrid along with a IS meson, for which the linear relations between matrix elements (103) and (102) are 
particularly simple: 

(({1I1P- ®lS) P ®L) lf \\V\\nl) = (104) 

(((inP + <8lS)p®.L)j / ||V||nO pa Jl{((lS®lP)p®L) lf \\S7\\nl). (105) 
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The expression (104) is a selection rule forbidding the decay of any £ state to a 1IIP~ hybrid in the limit that the 
hybrid quark radial wavefunction is a IP harmonic oscillator with the same size as the recoiling IS meson. The 
second expression contains a " «" sign to indicate that it is calculated with a first order expansion in the hybrid radial 
wavefunction. Further details of these approximations can be found in refs [4, 5]. 

In light of the experimental candidates 7Ti(1600), 7Ti(1400), the production of the exotic 1 *" along with a IS meson 
is an interesting possibility. A priori the most favourable production modes from ip and x states might be expected 
to have been 

ip -> 3 IUP- l 3 Si (106) 

Xo,i,2 -> 3 lIIPr 1 X S (107) 

In the proposed model the decay proceeds by the strong decay of a virtual state m qq in a S configuration. In this 
case the above modes should be strongly suppressed owing to a selection rule (104) above, as should the analogous 
production modes of the non-exotic J PC hybrids l 3 iIP^" 2 and ^EP^ belonging to the same family. Thus, for instance, 
observation of significant modes 

ifj -► TTlp (108) 

XI -* 7T17T (109) 

in the ir\ (1400), 7Ti (1600) channels would argue against a HIP - hybrid interpretation for them (see e.g. ref [4]). 

For charmonia states the non-relativistic approximation should be more robust, but there are questions as to 
whether the dynamics at this scale are driven by flux tube breaking. In general the selection rule presented here 
will be broken if the decay mechanism is dominated by perturbative gluons. On mass grounds it is possible that the 
X(3940) observed in 

e+e- -> V>X(3940), (110) 

contains a hybrid state with J h or J ++ quantum numbers. There has also been suggestion that the 1 F(4260) 
is a l-'-IIPT/ hybrid [7]: the immediate implication is that it should have (0,1,2) h partners l 3 nP^~j 2 , and notably 
the h and 1 h are expected to be lighter [10]. If e + e~ — > cc + cc can be modelled in a way analogous to the 
model presented here, a hybrid with J *" quantum numbers is not expected to be produced in the above reaction; 
however if the dominant production is via gluon exchange [3, 11] such hybrids could be produced. Predictions for 
such a mechanism require models that go beyond the present discussion. 

Subject to the proviso that production is by strong OZI, flux-tube breaking, rather than single gluon exchange, 
the selection rule may help discriminate between hybrid and other exotic or non-exotic interpretations of heavy or 
light-quark states, including those with non-exotic quantum numbers. 

The picture is altogether different for the production of the positive parity hybrids along with a IS meson. Con- 
centrating again on the hybrids with exotic J PC , the modes of interest are: 

^ -> 3 mp+ 1% (in) 

xo,i,2 - 3 inpJ 2 l 3 Sx (112) 



The relation (105) implies that to first order in the hybrid wavefunction there is a correlation of scale between the 
production rates of the above modes and those of conventional mesons. This relationship can be exploited to predict 
the branching fractions for the above modes relative to observed conventional meson modes. In general there is 
not an immediate relationship between the full decay amplitudes for any of the above 1IIP + IS modes and their 
conventional meson counterparts IS + ISP: charge conjugation requires that the modes have opposite spin so the 
angular momentum recouplings give different linear combinations of L, If matrix elements of the type (102) and (103) 
in the full amplitude. The ip decay mode above is a fortunate exception. It has been shown that the decays ip — > b\-K 
and ip ~ * a iP are consistent with the hypothesis that the virtual qq state through which the decay proceeds is a 3 Si 
state. For ip — > a,2p in particular the decomposition is uniquely in terms of the spatial matrix element on the right 
hand side of equation (105) (eqn. (82)), 

[ 3 SP 2 l 3 Si|n 3 S!] = i(((lP <g> 1S)i ® D)i||V||n 3 Si) 2 (113) 

while the analogous hybrid mode is expressed in terms of the spatial matrix clement on the left hand side of equation 
(105) [5] 

[ 3 mP+ l^oln'Si] = ^(((inP+ <g> 1S)i ® D) 1 ||V||n 3 Si) 2 (114) 

o 
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Using cqn.(105) there is a direct correlation of scale between the production amplitudes to first order in the hybrid 
wavefunction 

[ 3 mP+ l 1 S |n 3 Si] « i[ 3 SP 2 l 3 Si|n 3 S!] (115) 

not been made explicit here. 

In the following the implications for the isovector 2" 1 hybrid, denoted 62, will be considered; predictions for the 
isoscalar hybrid follow after adjustments due to flavour and phase space. The momentum of the b^ mode is very 
close to that of a-ip if the 62 mass around 2 GcV as expected, so that the above cqn. (115) translates into a direct 
relation between their branching fractions if the difference in the external phase space factors is disregarded. The a^p 
mode is second only to pn in magnitude and so even with the above suppression by a factor of 8 the corresponding 
hybrid mode 62"" should be observable: 

b.r.(ip -> 6 2 tt) ~ 1 x 1(T 3 (116) 

This result makes no reference to the initial state wavefunction other than that it is 3 Si. This first order estimate 
can only be taken as a guide as the leading order HIP wavefunction captures only the leading order angular 
dynamics. Corrections with a more realistic wavefunction are discussed in ref. [5] and it appears that the above is an 
underestimate. 

V. CONCLUSION 

Within the assumption that cc — > qq is the dominant intermediate state in the production of light hadrons, 
and that the strong OZI decay amplitude for qq —* qQ + Qq factorises in the sense of ref. [3], data imply that the 
intermediate qq state is not a hybrid meson. 

Although we have no well-developed dynamical model for this, it seems to us likely that this is consistent with 
the general expectation that hybrid qq production at short distances is suppressed, due to the nodes in both the 
qq and flux-tube wavefunctions. In turn this would imply that hybrid production from electromagnetic or weak 
currents, e.g. in e + e~ or B decays will be enfeebled. We note that the ^(4260), which has characteristics of hybrid 
charmonium[10, 12], has a nugatory leptonic width of O(cV), which was only exposed by study of the unusual channel 
■07T7T. If this state is not associated with hybrid charmonium, then even smaller leptonic widths would need to be 
accessed. Models or lattice QCD are needed to give insight into the short distance behaviour of hybrid wavefunctions 
and to assess whether they may couple dominantly through intermediate qQ + Qq loops. If the latter are important, 
this would go beyond our analysis, wihch has restricted itself to qq states. 

The production of hybrid via the long-range components of its wavefunction, such as OZI production in association 
with a conventional meson in cc decays, depends on the parity of the hybrid. Hybrids with negative parity are 
predicted to be suppressed; hence if either of iri (1400/1600) is a hybrid meson, we do not anticipate a significant 
signal in xi - ► ^^1 for example. Conversely, if cither is a qqqq member of a 10 or 10, there is no selection rule against 
their production. Hence a search for xi ~ * 7r7r i is merited; a strong signal would be interesting in its own right, as 
well as being a possible indicator that the exotic 1 *" ~K\ signal is not a hybrid meson. 

Our results suggest that positive parity hybrid production may be more promising. The prediction that -0 — » ~nbi has 
a branching ratio of ~ 10 -3 would make this typical in magnitude to other decays that have been studied successfully. 
The &2 is expected to occur in the 2 to 2.5 GeV region, where multibody decays may hinder its identification. However, 
the mode 62 — * 7ra2 is predicted to be a dominant quasi two-body channel [7] and hence this signal may be extracted 
in tp — * 57r. We would recommend a high statistics study of if) at BES in the hope of isolating this J PC exotic hybrid 
meson. 

This work is supported, in part, by grants from the Particle Physics and Astronomy Research Council, the Oxford 
University Clarendon Fund and the EU-TMR program "Eurodice", HPRN-CT-2002-00311. 
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